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Abstract. We investigate the relationship between ground-state (zero-temperature) quantum phase transitions in systems with 
variable Hamiltonian parameters and classical (temperature-driven) phase transitions in standard thermodynamics. An analogy 
is found between (i) phase-transitional distributions of the ground-state related branch points of quantum Hamiltonians in the 
complex parameter plane and (ii) distributions of zeros of classical partition functions in complex temperatures. Our approach 
properly describes the first- and second-order quantum phase transitions in the interacting boson model and can be generalized 
to finite temperatures. 



INTRODUCTION 

Quantum structural phase transitions (QPT's) at zero temperature are well known in both lattice [1] and many-body 
systems [2]-[24]. In any of such transitions, the structure of the ground state (and few low-lying states) does not evolve 
smoothly with the Hamiltonian control parameters, but flips abruptly from one configuration to another at a certain 
"critical point". A typical QPT Hamiltonian reads as 

H{X)=Ho + XV ={1-X)H{0) + XH{1) , (1) 

where Hq and V represent two incompatible terms, [IIo,V] ^ 0, and A is a dimensionless control parameter. Since 
arbitrary scaling and sign factors can be absorbed in V, one can require A G [0, 1]. The variation of A drives the system 
between two limiting A = and A = 1 modes of motions, which are in the many-body case typically associated with 
spherical and deformed shapes of nuclei [3, 7, 9, 11, 14, 15, 21, 23], paired and unpaired phases of strongly interacting 
Fermi systems [5, 6, 8, 17], or with normal and superradiant modes of some quantum optical systems [18]. 

For I A I large enough, the term AV in Eq. (1) prevails over Hq and the spectrum of //(A) just freely expands. The 
most interesting physics thus happens around the minimum Aq of the parabola that for finite Hilbert-space dimensions 
n determines the dispersion A^E — n^^trH^ — n^hr^H of energies in the spectrum as a function of A . The ground-state 
(g.s.) related structural phase transitions (if any) most typically appear at a critical value Ac not far away from Aq. They 
can be observed by analyzing the average 

{V)o^{Vo\vm = ^, (2) 

where £'o(A) is the energy and |^o(A)) the normalized eigenvector of the ground state. It can be shown that {V)q is a 
nonincreasing function of A, which in the QPT case develops either a discontinuity or nonanalyticity at Ac. The n^ °o 
limit is usually required for a QPT to occur (diagonalization of a finite-n Hamiltonian cannot generate a nonanalytic 
behavior), but distinctive QPT precursors can be often observed akeady in moderate dimensions. (Note, however, that 
a discontinuous change of the {V)o can take place even in finite-n cases if the ground state undergoes an unavoided 
crossing with another state [19].) If the (fc — l)th derivative of {V)q is discontinuous at Ac, the derivatives of Eq are 
discontinuous (singular) starting from the Kth one. This situation is described as a QPT of order K [2, 4], in analogy 
with the well-known Ehrenfest classification of thermodynamic phase transitions. 

Of particular interest are the situations when {V)o drops to zero and, simultaneously, also the dispersion, {V^)o, 
vanishes at the critical point. Typically, this may happen if V is semi-positively definite. Consequently, {V)q must 
remain zero for all A > Ac and the g.s. wave function gets fixed (in the non-degenerate case). In these cases {V)q may 



be considered as an "order parameter" that distinguishes two quantum "phases" of the model (the values {V)q = and 
(y)o 7^ being attributed to "more symmetric" and "less symmetric" phases, respectively). 

The ground-state QPT's happen at temperature T = and thus have no real thermal attributes. Therefore, questions 
often arise whether the term "phase transition" used in this context is just a metaphor, or whether it refers to some 
deeper analogies in standard thermodynamics [25]. Of course, the situation becomes more general if the system is 
described with both A and T taken as free parameters. The g.s. average {V)o is then replaced by a thermal average 

_ tr(y-.-^-'^) dFo 
^^^^= tre-T-^H ^^T' 

■where Fq{?c,T) = —Tin tve^^ ''^f'''' is the equilibrium value of the free energy and V' = —r[^e^^ '^(-^jje^ 'hW 
an operator generating the finite-temperature "order parameter" (one can prove that {V')t {V)o for T 0). In this 
case, both structural (driven by A) and thermodynamic (driven by 7") phase transitions turn out to be just two different 
aspects of the same phenomenon, residing in the full A x T parameter space [1, 2], and a common language should be 
developed for their description. 

In this contribution, we will mostly deal only with the simplest case, i.e., with the T = limit of quantum phase 
transitions. The approach [24] will be presented that makes it possible to treat such transitions fuUy in parallel with 
thermodynamic phase transitions at finite temperatures. This approach is based on a similarity between the distribution 
of zeros of the partition function Z at complex temperatures for systems undergoing classical phase transitions [26, 27] 
and the distribution of so-called branch points of QPT Hamiltonians (1) in complex-extended A plane [6]. Since a 
generalization of our method to finite temperatures seems possible, we believe that it represents the right way toward 
the unified description of quantum phase transitions in the A x T space. 



QPT'S IN THE INTERACTING BOSON IMODEL 

As a testing ground, we will use the interacting boson model (IBM) [28], well known from nuclear- structure studies 
as well as from analyses focused on general features of quantum chaos and quantum phase transitions. This model 
describes nuclear shapes and collective motions in terms of an ensemble of A' interacting s and d bosons with angular 
momenta and 2, respectively. Its algebraic formulation is based on the dynamical algebra U(6) formed by 36 bilinear 
products b^bj of boson creation and annihilation operators, where bj = s or df^ with fi = —2, . . . , +2. This allows one 
to extract several alternative dynamical symmetries and find analytic solutions for the corresponding Hamiltonians. 
All the IBM dynamical symmetries lead to the algebra 0(3) of angular momentum L = ^/lO{d'^d)^^^ which guarantees 
the invariant symmetry of the model under rotations. 

Individual dynamical symmetries are named after the first algebras in the respective U(6)D . . . dO(3) decomposi- 
tions, so we have U(5), SU(3), 0(6), SU(3), and 0(6) dynamical symmetries (the fourth and fifth chain differs from 
the second and third one, respectively, just by relative phases between s and d bosons in definitions of the correspond- 
ing algebras). Both first- and second-order QPT's — according to the Ehrenfest classification — are present in the IBM 
parameter space between these symmetries. Since geometry can be attributed to given algebraic structures via the 
coherent state formaUsm, the IBM phase transitions describe changes of the nuclear ground-state shapes. In particular, 
using the projected coherent states [29] 

|iv,^,r)- (/+/3cosr4 + ^[^l2+42])''|o) (4) 

with j3 and y interpreted as Bohr geometric variables, one obtains the following shape types associated with individual 
dynamical symmetries: spherical [U(5)], prolate [SU(3)], oblate [SU(3)], and deformed y-soft [0(6), 0(6)]. 

The IBM Hamiltonian is assumed to have a general form containing one- and two-body terms, with few restric- 
tions resulting from fundamental symmetry requirements. This leads to the most general Hamiltonian with six free 
parameters (except an additive constant). However, a more involved analysis focused on phase-transitional proper- 
ties [7] reveals that there are only two essential parameters. An archetypal two-parameter Hamiltonian (see, e.g., 
Refs. [10, 13]) reads as 

QxQx' 



H^{X) = {\-X) 



N 



-And, (5) 



where the parameters change within domains A e [0, 1] and % G [-3^, ^y^], while nd=d^ -d represents the d-boson 
number operator and Q^^ = d^s + s^d + xid'^d)^^^ the quadrupole operator. The U(5) dynamical symmetry is located 
at 77 = 1 and X arbitrary, 0(6) at {ri,x) = (0,0), SU(3) and SU(3) at {ri,x) = (0,-^) and (0, + ^), respectively, 
while 0(6) is absent in the given parametrization. 

Clearly, if x is fixed to a constant, the Hamiltonian in Eq. (5) has the general form (1) with 

y=nd + ^^^ (6) 

being a semi-positive operator. It is not difficult to see that {V)o > for X ~0 and {V)o ~ for A = 1, so the scenario 
described in the previous section may be expected to apply in the N ^ 0° hmit. Indeed, the order parameter can be 
expressed in terms of the Bohr deformation parameters JSq and resulting from minimization of the coherent- state 
energy functional 

1 (5A -4)/32 +4^1(1 - A)/33cos3r+ [A - '^xH^ - A)]/?^ 

(normaUzation per boson is included to ensure finite AT — > 00 asymptotics), namely 

For ;^ ^ 0, the value of j3o changes from /3o 7^ to at 



4 + fz2 



^ = ^-T^ > (9) 

indicating a first-order deformed-to-spherical QPT (the order parameter (l^)o drops from a positive value to zero). 
For ;|f = 0, the value j3o °^ VAc — A valid in the left vicinity of Ac continuously joins with j3o = vahd above A^; 
the corresponding QPT is of the second order, with the critical exponent for (l^)o equal to 1. Because of the proper 
weighting of both terms in Hamiltonian (5) the finite-A^ precursors of the QPT behavior are always located in the 
region around (9), independently of A'. 

Besides the spherical-deformed phase separatrix, there exists also a separatrix at x = and A < | which corresponds 
to the first-order QPT between prolate and oblate shapes, where either j3o changes the sign or, equivalently, 70 jumps 
from to ^. In this case, however, the form (1) [with A replaced by x, and Hq and V by the corresponding expressions 
following from Eq. (5)] can be used only locally, close to ^ = 0. We will not be dealing with the prolate-oblate 
transition here. 

An obvious way to classify the IBM deformed-spherical QPT (driven by the interaction parameter A) relies on the 
analogy between theA' ^ °o g.s. energy per boson, S'o{X,x) = <^(A,x;j3o)?b)> as a function of A, and the equihbrium 
value of the free energy, Fo{T), as a function of temperature T. This leads to the Ehrenfest classification. Using the 

standard definition C = —T-^Fq, one can even introduce a QPT analog of the "specific heat" 



^ = -A^ = -A^^= lim^;^y (10) 

and show that it behaves exactly as expected for a thermodynamic phase transition of the respective order [16]. Here 
and in the following, = Ei{X,X) stands for the ;th eigenvalue of H^iX) while {V)ij = {^i{X,x)\V\^j{X,x)) and 
(y), = (y)„ denote matrix elements of y involving ith and yth eigenstates. Note that the suminEq. (10) runs effectively 
only over those states that have the same synnmetry quantum numbers as the ground state, since otherwise {V)oi = 0. 



ALTERNATIVE DEFINITIONS OF A QPT "SPECIFIC HEAT" 

It turns out that Eq. (10) does not represent the only form of the QPT analog of specific heat. Inspired by the thermo- 
dynamic relation C = T-^S, where S is the entropy, one can define = A with = - 1:,- 1 {i\Wo) pin | {i\Wo) p 
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FIGURE 1. "Specific heat" (13) with n = 1 including all 7 = states for (a) the first-order and (b) second-order QPT of the 

interacting boson model [Hamiltonian (5) with (sl) x = and (b) X = The curves in order from the lowest to the highest 
correspond toN = 10, 20, 40, and 80. The insets show the increase of the peak maximal amplitude with A'. 



being the wave-function entropy of the ground state with respect to the eigenbasis |/) of //(I), i.e., the basis associated 
with the U(5) dynamical symmetry in our case. This seems to be a plausible alternative definition of specific heat [16] 
for systems where {V)o and (y'^)o drop to zero at Ac. (Note that A should be inverted here to A = 1 — A in order that 
the entropy y increases with "temperature" A.) 

Another possibility is to randomize Hamiltonian (1) by adding a small stochastic component 5A to the control 
parameter (with (5A) = and (SA'^) = cr^ <C 1) and to exploit the perturbative expansion H{X + 5A) = H{X) + 
(5A)y. In this way, the ground state wave-function transforms into a density operator Po(A), whose identification 
with a canonical ensemble provides a natural basis for a noise-induced thermalization of the given quantum system. 
The resulting specific heat reads as 



' = tr(poln2po) - tr2(polnpo) « [a^ln^ a^) £ 

i>0 



\{V} 



•Oil 



(11) 



where the last equality is valid for very small values of a (in the IBM, a natural scaling of a for ^ oo is 
oc [16]). Again, formula (11) evaluated for the IBM Hamiltonian (5) is peaked just around Ac, which indicates 
an increased mixing of the Hamiltonian eigenfunctions in the QPT region. See Refs. [12, 16] for details. 



The last analog of specific heat will be elaborated in the rest of this contribution. It originates in the expression 



-£ln|£i-£o| 

i>0 



(12) 



for the "potential energy" of a given set of levels with the same symmetry quantum numbers, as derived from 
the relation between one-dimensional distributions of charges in a planar universe and analogous distributions of 
eigenvalues in Gaussian matrix ensembles [30|. This relation, referred to as the static Coulomb-gas analogy, must 
be distinguished from the dynamical Coulomb-gas analogy, introduced by Pechukas and Yukawa [31], that describes 
eigenvalue dynamics for Hamiltonians of the form (1). 

A quantity proportional to ^ in Eq. (12) can be considered as another QPT analog of the thermodynamic potential. 
In particular, if is associated with = ii~'A^, where i2 is a scaUng constant that is to be discussed later, the 
specific heat takes the form 



2A^ 
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with the last expression resulting from the Pechukas- Yukawa equations [the last two expressions in Eq. (13) are given 
here just for comparison with Eqs. (10) and (1 1); from the computational viewpoint the most convenient definition is 
of course represented by the first formula]. Since the motivation for the proposed relation between and ^ may 
be unclear at this moment, we ask the reader for patience till the next section, where the foundation of Eq. (13) will 
become clear. 

It turns out [24] that the specific heat given by the last formula exhibits a form that is very similar to that of Eq. (10). 
This is illustrated in Figure 1, where specific heat (13) with i2 = 1 is shown for the IBM first- and second-order phase 
transitions. Since the AT — » oo limit is numerically inaccessible, we present here only calculations for A'^ = 10, 20, 40, 
and 80, corresponding to the curves with increasing height, the maximal values of '^i' being shown separately in the 
log-log insets as a function of A'. Clearly, the peaks in panel (a) are sharper and higher than those in panel (b), as 
indeed expected in the first- and second-order phase transition. Moreover, as will be argued below, the inclusion of the 
right dependence of the scaling Q. on dimension n leads to the correct A' ^ °o asymptotics of Eq. (13). In this case, the 
specific-heat values at A^ converge or diverge, respectively, for the IBM second- or first-order QPT's. 

All the above-proposed definitions of specific heat lead to certain sums of squared interaction matrix elements 
divided by some combinations of energy differences. The observed peaked behavior of these quantities results from 
the fact that the spectrum becomes more compressed and, simultaneously, also the sum of level interactions get stronger 
in the phase-transitional region. In fact, we are dealing here with a phenomenon of multiple avoided crossing of levels 
[10]. The analytic description of such effects can be given in terms of so-called branch points of Hamiltonian (1) in 
the complex A plane [32]-[36]. As will be shown below, this description is particularly relevant for the last definition 
of specific heat in Eq. (13). 



BRANCH POINTS AND ZEROS OF PARTITION FUNCTION 

Branch, or exceptional points are places in the complex plane of parameter A where various pairs of eigenvalues of 
the complex-extended Hamiltonian (1) coalesce [32]. They are simultaneous solutions of equations dst[E — H{X)] = 
and det[E — H{X)] = 0, that after elimination yield the condition [33, 35] 

^ = U^k = H-^U(Ej-Eif = (14) 

k i<j 

with = Tli{ytk){Ei — Ek). The discriminant ^ is a polynomial of order n{n — 1) in A (the dimension of the Hilbert 

space n being now assumed finite) with real coefficients and its roots thus occur as complex conjugate pairs. 

Except at these points, the complex eigenvalue E(X) obtained from the characteristic polynomial of Hamiltonian (1) 
is a single analytic function defined on n Riemann sheets. The energy labels in Eq. (14) enumerate the respective 
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FIGURE 2. Branch points corresponding to all 7 = states for (a) the first-order (|;|;| — ^) and (b) second-order {% — 0) QPT 
paths of the IBM Hamiltonian (5) with A' = 20. Diamonds represent single branch points and crosses the degenerate ones (within 
given numerical precision). 



Riemann sheet according to the ordering of energies at real A. The degeneracy points are square-root branch points 
where the Riemann sheets are pairwise (in generic cases) connected. The leading-order behavior on the two connected 
sheets close to the branch point Ab is given by E{X) —E{Xh) ~ a^/X — Aq (as a doubly-valued function), with a being 
a complex constant [6, 34, 35]. 

The relation of branch points to QPT's has been declared several times — see, e.g., Refs. [6, 10]. Clearly, a branch 
point located close to the real A axis affects the local evolution of the corresponding pair of real energies so that the 
two levels undergo an avoided crossing with accompanying rapid changes of wave functions. A cumulation of branch 
points close to some real point A^ thus can give rise to massive structural changes of eigenstates, as observed in QPT's. 
It moreover turns out that the density of branch points close to A^ may be determinative for the QPT order, as observed 
on the real A axis. 

The simplest example can be found in dimension n = 2. A general real 2x2 Hamiltonian of the form (1) reads as 

It is easy to show that the spacing between both eigenvalues Ei — = AE behaves according to \AE\ = 
VX+BAT+CA?, where the coefficients A > 0, B, and C > are some combinations of constants in Eq. (15) 

such that the value D = — 4AC is negative. The minimal distance of both levels, \AE | = Aq = \J is reached at 

A=Ao = — ^.IfAis extended to the complex plane, one finds that both complex eigenvalues Ei and £2 cross in two 
complex-conjugate branch points Ab± = Ao ± ; Aq. For Aq converging to zero, the branch points come to the real axis 
and the avoided crossing turns into a real crossing, implying lA^I = \/C(A — Ab). The line that passes both points Ab+ 



and Ab- in the direction perpendicular to the real axis represents the locus of points where either real or imaginary part 
of El — E2 vanishes. In particular, one has ReAZi = for |lmA | > |ImAb± | and ImAE = for |ImA [ < |ImAb± [. Only 
at A = Ab± both loci overlap so that at these points ReA£ = ImA£ = 0. This corresponds to the above-mentioned 
behavior of the complex square root (vaUd for Aq > 0). 

In the n = 2 case, the domain of E{X) consists of two Riemann sheets only. It is evident that the analysis becomes 
much more complicated for higher dimensions. Methods to locate individual branch points in the complex plane for 
arbitrary n were developed (see, e.g., Ref. [33]), but without the knowledge of explicit Riemann sheet structure of 
the specific problem one cannot select those branch points that are located on the ground-state sheet. An example 
is shown in Figure 2. Here, branch points corresponding to the 7 = submatrix of the IBM Hamiltonian (5) with 
A'^ = 20 are shown for the first-order QPT with % = ^^ (panel a) and for second-order transition, X = ^ (panel b). 
The distributions in both panels show no cumulation of branch points near the respective critical points Ac = 0.8181 . . . 
(a) and = 0.8 (b), which reflects the fact that all Riemann sheets are mixed together. 

As can be observed, some of the branch points in panel (a) of Fig. 2 come very close to the real axis, but the 
absence of additional integrals of motions in the X 7^ transitional regimes guarantees that ImAb 7^ in all these 
cases. On the other hand, several branch points in panel (b) are located exactly on the real axis, with ImAb ~ 0, which 
seems to contradict the familiar no-crossing rule. Indeed, the % — ^ transitional path is not quite generic since the 
underlying 0(5) dynamical symmetry of the 0(6) and U(5) hmits remains unbroken and gives rise to the seniority 
quantum number v, that is valid along the whole transition [37]. Consequently, levels with different v's can cross 
without repulsion. It must be stressed, however, that none of these very close avoided (panel a) or unavoided (panel b) 
crossings is related to the ground-state Riemann sheet. 

In generic cases, the determination of the Riemann sheet structure in its full complexity is prohibitively difficult 
even for moderate dimensions. The loci of ReAZ? = and ImzVE = crossings are not any more straight lines, as 
for n = 2, but generate comphcated patterns that can be hardly disentangled having only a finite numerical precision. 
(Nevertheless, the ReAE = curves for IBM form a flow with prevailing perpendicular orientation toward the real A 
axis, see the preprint in Ref. [24].) As a consequence, practically nothing is known about the density of the ground-state 
branch points in a vicinity of A^ for QPT's of various orders. 

We will show that the specific heat in Eq. (13), although depending solely on the real-A observables, represents an 
indirect measure of the density of branch points on the ground-state Riemann sheet near the real axis. To prove this, 
we start from following question: If ^0 = ii^'A'^ represents the equilibrium value of the "free energy", as assumed 
above, what is the corresponding "partition function" 2f? Using the thermodynamic relation Fq = —TlnZ, one finds 

^" = n(S,-Eo), (16) 

!>0 

thus ^ coincides with the H^'th power of the partial discriminant ^0 from Eq. (14). Recall that the square S'^ is 
a polynomial in A with n—l complex conjugate pairs of roots, each of them being simultaneously assigned to one 
other i^^, . These roots correspond to the branch points located on the A:-th Riemann sheet. Thus branch points on the 
ground-state Riemann sheet are zeros of the fictitious partition function iF, that generates the specific heat in Eq. (13). 

Zeros of canonical or grand canonical partition functions Z in complex temperatures and/or chemical potentials play 
an essential role in the fundamental theory of thermodynamic phase transitions. Indeed, as noticed by Yang and Lee 
[26], for a system of particles with mass m interacting through a potential with a hard repulsive core the grand partition 
function Z can be written as a polynomial in the parameter y — (/i^ ' VmT)^ exp (T^ '/i), containing both temperature 
T and chemical potential 11. While for a finite size of the system the roots of Z{y) must keep away from the real y axis, 
in the thermodynamic Umit there may exist places where the zeros approach infinitely close to the real axis, giving rise 
to a phase transition at the given value jg. This method was adapted also for canonical ensembles [27], where it was 
shown that the density of zeros of the partition function close to the phase-transitional temperature Jc even determines 
the order of the transition, as reflected by the behavior of standard thermodynamic specific heat. 

This was recently proposed [38] as a basis for phase-transitional analyses in small systems. For instance, if zeros 
of Z are located along a curve crossing the real axis and if the closest zero converges to a real point in the 
thermodynamic hmit, the order of the corresponding phase transition at is determined by (i) the power a in the 
dependence p (ImF)" of the density of zeros close to the real axis, and by (ii) the angle v between the Z = 
domain and the normal to real axis [38]. Namely, the phase transition will be of the first order for a = V = 0, of the 
second order for a € (0, 1], and of a higher order for a > 1. A detailed analysis of even more possibiUties can be found 
in Ref. [27]. 

Because the "specific heat" (13) results — as we know now — from the "partition function" (16), it basically measures 
the density of zeros of 3f near the real A axis, or, in other words, the density of branch points on the ground-state 



Riemann sheet. The peaked behaviors shown in Fig. 1 thus indicate that the IBM g.s. branch points indeed accumulate 
near Ac for N ^ °°, and that the degree of this accumulation differs for transitions of different orders. We are now 
in a position to formulate the central statement of this contribution, namely, the surmise [24] that the Kth-order QPT 
distribution of the ground-state branch points is quantitatively similar to a distribution of the Z{T) complex zeros in a 
thermodynamic phase transition of the same order K. 



NORMALIZATION AND LARGE-N ASYMPTOTICS 

To test the conjecture proposed at the end of the last section, we need to discuss the factor Q,, that appears in Eqs. (13) 
and (16). In fact, this factor should not be a constant, but a certain function depending on the system's size (quite 
similarly to the standard specific heat, which must be normalized to a unit part of the thermodynamic system). Since 
the polynomial is determined (up to a multipUcative constant) by its complex roots, one may express specific heat 
(13) as an integral containing the density Po(A) = 5(A — Ab,) of the branch points Ab, on the g.s. Riemann sheet 
[24] . The normalization of / podX to unity yields a factor o= (n — 1 ) ^ ' (there are n—1 branch points on each Riemann 
sheet) and one arrives at the formula 

n = no(«-i), (17) 

where n is the relevant Hilbert space dimension and an arbitrary constant. Eq. (17) is analogous to the scaling based 
on the grand partition function, as discussed by Yang and Lee [26], since in that case the order of the polynomial Z(y) 
equals to A^max, the maximal number of particles in a given volume, and Q. oc Nm^x can be identified with the volume. 

In the IBM, the total number of states grows roughly as ~ A^^ / 120 for very large boson numbers (if not counting the 
rotational degeneracy). However, the dimension of the / = subspace, which is relevant in the calculation leading to 
Fig. 1, is given by n ^ N^/12. So the correct normalization of Eq. (13) in the IBM case is by a factor Q.°^ N^. Figure 3 
shows — by the curve demarcated by squares — the maximal values of such normalized specific heat in the second-order 
(X = 0) QPT region for very high boson numbers A'^. Note that such high-dimensional calculation were enabled by 
the underlying 0(5) dynamical symmetry at j = 0, so they could not be performed for the other transitional path (a). 
The algebraic increase shown for < 80 in the inset of Fig. 1(b) (calculated for O = 1) represents only the initial part 
of the curve in Fig. 3, where the convergence to the A oc asymptotics is evident for A'^ > 300. For comparison, the 
convergence of the maximal value of specific heat (10) to its asymptotic value (which is for X = equal to 12.5) is 
also shown in Fig. 3 by full circles. Clearly, the degree of convergence in both cases (squares and circles) is about the 
same. On the other hand, since the A^ < 80 increase observed in panel (a) of Fig. 1 is much faster than the increase 
in panel (b), one expects that the asymptotic increase of the maximal value in the first-order QPT is faster than 
Ao< (although in this case the A'^ > 100 region is numerically inaccessible). This indicates that the QPT behavior 
of the specific heat (13) is consistent with the behavior of the standard specific heat in the first- and second-order 
thermodynamic phase transitions, in agreement with the above proposed conjecture [24]. 

This conclusion is also supported by the analysis of the A'^ — > «> Umit. Indeed, due to the degeneracy of j3o = 
and j3o 7^ minima of the function (7) at Ac for J 7^ 0, there will be a branch point on g.s. Riemann sheet that 
asymptotically reaches the real axis. In contrast, there is no actual ground-state involving degeneracy in the A' ^ °o 
limit of the second-order transition, although, as can be shown, exceptional points come infinitely close to the real axis 
at the critical point. 

Let us stress that in calculations leading to the above-discussed results all levels with / = were included in Eq. (13). 
We know, however, that for X = some pairs of levels actually cross at real A (due to the seniority quantum number), 
which implies that first and second derivatives of the corresponding energies may be discontinuous or diverging. 
Fortunately, a detailed consideration of this problem reveals that all these singularities cancel out exactly and do not 
affect the shapes of curves in Figs. 1(b) and 3. Indeed, if £, crosses with £,+1, the discontinuity (singularity) of the first 
(second) derivatives in Eq. (13) for both levels is the same, but with opposite signs, so the sum of both contributions is 
incremented correctly, as if the levels continuously passed the crossing. 

In any case, one can repeat all j = calculations including only the v = subspace of states with / = into 
expression (13). Since this subspace contains the ground state and does not mix with v j^Q subspaces, one obtains a 
plausible new definition of the X=^ specific heat, which avoids the above-discussed problem with singularities. The 
dimension of the v = subspace grows as n ^ A'/2, thus i2 «; A^. Maximal values of the resulting specific heat curves 
with this normalization are shown in Fig. 3 by triangles. We observe a similar dependence as before, but because of 
lower dimensions the convergence with A'^ is much slower now than before, not allowing a fuUy conclusive evidence 
of finite asymptotics. The available data for boson numbers up to 5000 (see the inset in Fig. 3) are consistent with 
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FIGURE 3. Maximal values of specific heat (13) (squares and triangles) normalized according to Eq. (17) in the IBM second- 
order QPT (x = 0) for very large boson numbers. Squares (or triangles) correspond to the inclusion of the whole 7 = space (or its 
V = subspace) into the sum in Eq. (13). Different (arbitrary) scaling factors were used to show both sets of data within the same 
range. For comparison, the convergence of the maximal value of specific heat (10) to its asymptotic limit = 12.5 is shown by 
full circles. 



large-A'^ behaviors given by A oc A^* for k G (0, 1.35), in agreement with the expected value k=l. Let us stress that 
the dimension n of the V = subspace atN = 5000 is comparable with the dimension of the whole ./ = subspace at 
A' ~ 170, so the absence of saturation in the v = case in the available domain of boson numbers is not surprising. To 
obtain v = results comparable with 7 = at A/^ ~ 300, one would need to go up to A/^ « 15000, which is beyond our 
present computational possibilities. 



CONCLUSIONS AND OUTLOOK 

We have discussed several aspects of quantum structural phase transitions at zero temperature that could eventually 
lead to their unified description with classical thermodynamic phase transitions. The most fundamental aspect seems to 
follow from the analogy [24] between the branch points related to the ground-state Riemann-sheet of Hamiltonian (1) 
and complex zeros of the partition function of a classical phase-transitional system. Numerical data strongly supporting 
this conjecture were obtained in the interacting boson model, although a fully conclusive test would require to further 
increase the upper limit of available boson numbers. The immediate next task is to extend the calculations presented 
here to other systems that exhibit quantum phase transitions of various orders. 

An open problem is the generalization of our approach to finite temperatures. In this case, both structural and 
thermodynamic phase transitions describe nonanalytic behaviors of the thermodynamic potential Fo{X,T) along two 
perpendicular directions in the XxT parameter space. While the T-direction is described by standard thermodynamics, 
involving the theory of complex zeros of the partition function, the A-direction is not as familiar. Formula (3) represents 
a possible starting point of the analysis. The finite-temperature average of V is just a weighted sum of averages 
corresponding to individual excited states 

{V')T = ^Z^-'"'^-{V')i (18) 

^ i 

and one can ask whether a thermally weighted sum pr of branch-point densities p, on the Riemann sheets assigned to 
excited states will, in the finite-r case, take the role of the g.s. density po, as discussed above. 

The interacting boson model seems to be an ideal tool for studying various aspects of complex quantum dynamics, 
including quantum phase transitions. One of its advantages in this respect is the presence of quantum phase transitions 



of both first and second orders. We consider this model as one of the best candidates also for future finite-temperature 
studies. 
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